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The radiation fields of a magnetic line source axially oriented on a conducting cylinder 
which is in a plasma environment have been analytically determined. The plasma environ- 
ment is taken to be a series of homogeneous cylindrical layers none of which is contiguous 
with the cylinder. It is demonstrated that a quite significant distortion of the field pattern 
is effected by the presence of a lossless plasma sheath. Characteristics of the radiation 
fields are discussed and radiation patterns, computed in a plane through the axis of the 
cylinder and the line source, are presented. The patterns are found to have large narrow 
amplitude fluctuations that can be attributed to leaky waves. 

1. Introduction 

The gross effects on radio wave propagation of the plasma sheath which envelopes hyper- 
sonic vehicles are fairly well understood. It is recognized that when such vehicles travel at 
sufficient speeds in the upper atmosphere there may be a blackout of radio communication and 
sufficient distortion of radar signals and patterns to make radar systems unreliable. 

The most direct approach to increasing communication and radar system reliability when 
the sheath environment is present is to increase the propagation frequency well beyond the 
peak plsama frequency in the shock wave. In this way the sheath is made to appear trans- 
parent to the waves so that, aside from an increase in the probability of breakdown on transmit, 
the sheath might just as well not be there. The problem is, however, that the plasma frequency 
in some regions of the sheath can go well beyond the frequency of operation of present systems. 
In regions of the sheath aft of the nose of a vehicle there is a significant decrease in electron 
density so that, in general, there is a point at which an rf system will operate with only modest 
distortion of the propagating signals. It is such a region that is of concern in this paper. In 
particular, we shall investigate the effect of the sheath on propagation when the propagation 
frequency is greater than the plasma frequency, but when they are both of the same order of 
magnitude. 

The geometrical configuration of the hypersonic plasma sheath in terms of the spatial 
variation of both charged and uncharged particles has been subjected to a great deal of investi- 
gation over the past few years, and while an exact determination of the spatial variation of the 
constitutive parameters of the sheath has not been made for even the simplest of vehicle con- 
figurations, a general description is available for cylindrically symmetrical vehicles with tapered 
nose sections. In this report the radiation characteristics of a linear slot antenna flush mounted 
to a circularly cylindrical vehicle will be considered. In the sections 2, 3, and 4 the enveloping 
plasma sheath will be taken to be a series of consecutive plasma layers contiguous with each 
other but separated from the body of the vehicle. This separation is characteristic of hyper- 
sonic shock waves at other than stagnation regions. The vehicle, as well as the shock wave 
region, will be taken to be infinitely long and the cylinder will be considered to be infinitely 
conducting. In section 5 an approximation technique for determining the field distortion due 
to sheaths which are not necessarily cylindrical and composed of homogeneous layers will be 
developed and discussed. It will be shown that this technique predicts the significant distor- 
tions which have been computed for the single layer case under the assumption of leaky waves 
excited by the source. 



iThis work was sponsored by the Electronics Ileseaich Directorate, Air Force Cambridge Research Laboratories, under Contract No. 
A F 19604-8386. 
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2. Formulation of the Radiation Problem 

2.1. Solution for the Field in Cylindrical Regions 

The radiation field of a magnetic line source lying on a circularly cylindrical conducting 
surface is to be determined. The surface is infinite in extent, but the line source is finite and 
lies parallel to the cylinder axis. Encircling the surface is a series of concentric cylindrical 
regions containing homogeneous plasma media with arbitrary electron densities, except for the 
inner and outer regions which have zero electron density. This geometry is depicted in figure 1 . 

The technique for obtaining a solution to this cylindrical radiation problem is well known. 
The solution may be found in terms of E and H modes referred to the z direction, which forms 
the cylinder axis. Expressed in terms of potential functions, F and A, satisfying the scalar 
Helm hoi tz equation the two sets of modes 



E 1 = -i W /iVXW H 2 =icoeVX(^U 2 ) 

are defined for each cylindrical region. The total field is then 



(i) 



E=E,+-!- VXH 2 H=H 2 -^- VXE, 
^coe ^co/x 



(2) 



The propagation and plasma frequency are considered to be in the kilomegacycle region, 
so that the cyclotron frequency, which is entirely due to the earth's magnetic field in the absence 
of sources of low frequency magnetic field on the vehicle, may be neglected. The constitutive 
parameters in the cylindrical regions are thus isotropic. Further it is assumed that each region 
is characterized completely by an equivalent constant permittivity. 

Solutions for A and F may be given in terms of cylindrical wave functions in the form 
[Wait, 1959a]. 



A= S 



F= S 



f- 



dh{a m (h)H%Uk 2 -h 2 P )+A m (h)J m (JW^h 2 p) } ■ 



-i(m<t>+hz) 



dh { b m (h)H%\^[W=h i p)+B m (h)J m (W = Vp) } «-«■♦+*> 



(3) 



and the fields may be found directly from eq (2). The form for the field components transverse 




Figure 1. Geometrical configuration for cylindrical 
plasma sheath problem. 
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■iim^-^hz) 



to the radial direction in any region is 

m= — co J -a, 

H z = 2 f" rf/^(F-/r)[6„,(/0^ 2) (VF=r 2 p)+ J B„,(A)J m (VF=^p)]}e- i< '"*+ ftz) 

m=—CQ J -CD 

m = -co J -co v^ 

r°° r «p , , . 

TO =-co J_oo ^ M^ 



mh , 



'}• 



^^[b m (h)H^(^k 2 -h 2 p)+B m (h)J m (^k 2 -h 2 P )} Ye-'***". (4) 
P 

Equations (4) represent the complete solution for the field in any homogeneous region of 
space with characteristic propagation constant k. The set of coefficients a mj A m , b rny and B m 
are determined for each region by matching the four field expressions at eacli boundary. Since 
the fields are finite in each region tin 4 transform operation may be considered to be a factorable 
operator and the matching technique consists of equating the quantities inside the braces 
in eq (4). 

2.2. General Formulation, for the Multilayer Problem 

To facilitate the writing of solutions the convention adopted will be that the cylindrical 
surfaces are at radii pj where the cylinder surface is p and the outer plasma region is hounded on 
the outside by p N . The regions themselves are numbered the same as their outer boundaries 
so that /^represents the propagation constant of the /th region where / has t he values 1,2,. . . , 
N-\- 1 and k is reserved for the plasma free regions, i.e., 

The expansion coefficients of the fields are labeled with a superscript appropriate to their 
regions in the form a$ and as in [Wait, 1959a] a new symbol is defined: 

Uj=^k 2 j—h 2 . 
Equating the braces of (4) at the/th boundary yields the four equations 

-^ [a^H^(n J+lPj )+A^J m (u j+lPj )] 

pj 

jUCO Pj 



ih 2 



fXiO 
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An analytic solution to the problem is found by determining the 4(iV+ 1) modal coefficients 
a\l\ A%>, b%>, BH ) ;j=l, 2, . . ., N+l. Equations (6), when written for the N surfaces, repre- 
sent 4iV equations. The necessary additional four equations are obtained from the radiation 
condition which requires that there be only outgoing waves in the external region and from 
matching the fields to the source at the cylinder surface. Thus, 

ir^^r^o (7) 

and, at p=p 

E z = S P dh{u\[a^H m (u 1 p a )+^J m (u lPo )]}e- t ^+'">=0 (8) 

171 = — co J -co 



E*= 2 f" dh\i mUl [b£H^'(u iPo )+B^J' m (uiP<,)] 

m= — co j _oo I 



CO (* oc 

TO=-co J - 



Po J 

dh{P m (h)}e- i(m * +h2) (9) 



where the expression on the right side of (9) represents the expansion of the line source in a 
Fourier integral in z and a Fourier series in </>. Equation (9) may then be simplified by equating 
only the braces. 

2.3. Asymptotic Evaluation of the Radiated Field 

To find the radiated field it is necessary to carry out the integration and summation indi- 
cated in (4) when < /=A7'+l. The usual technique is to change the cylindrical field coordinates 
p, <t>, and z to the spherical coordinates R, 0, </> (see fig. 1) and to evaluate the integral asymp- 
totically for large R. Thus under the substitutions 

p=R sin 

z=R cos 

the <t> component of the electric field is 



CO /*oo f 

#♦= S dh-l iuixu N+1 b^ +1) H^'(u N+ ,R sin d) 

m = — co J _oo v. 



a ( m N+ vH%\u N+1 R sin 0) \ 6 -«"*+»* «>*<»• (10) 



R sin 



For finite m and large value of the argument Uv+iR sin 0, the Hankel function is approximately 
[Wait, 1959b] 



ff£> (u N+1 R sin 0) « J ^— e-t k+i* sln '-^ *] 

\ wu N+1 R sin 



(id 



' irU N +iR sin 
and its derivative is 

H$'{u N+l R sin d)=H^ 1 (u N+1 R sin 0) ^— - H^(u N+1 R sin 0) 

Uj^j^iiX Sill C7 

V2 r 2m + 1 1 

„ . „ 6-*K+i* Blne — 4" »]. (12) 
ttu n+ iR sin 

Substitution of (11) and (12) into (10) yields 

y^J / ? -ifw^+jiSsinfl-^ — tt] I -i^+hR cos 0> Qg\ 

\ \ tu n+1 R sin J 
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which is an integration on the real axis of the A-plane. The radiation field is the sum of terms 
of order I JR. Terms of this order are readily determined from (13) by evaluation of the integral 
by the method of stationary phase which yields as solution for the integral [Copson, 1946] 



the quantities 



i=i 



1= r ' e iR t> M 4<(x, R)dx 
J* -« 



(14) 



where the upper solution is obtained if x is not a stationary point and the lower solution holds 
when x is a stationary point. Examination of (13) clearly indicates that the term for which 
a ( * +1) is a factor yields no radiation fields while the b ( J? +1) factor yields a radiation field only 
at the stationary points of the phase term 



which are at 



(h cos d-\-u N+l sin 6) = (h cos d+^kl—h 2 sin 0) 
h=k cos 6. 
From (16), (15), and (13), then, the radiation field is given by 



E $=21^11 



s,—ikR a> . /t A 



R 



S «*+i&2r i) We' 



h = kQ cos i 



(15) 
(16) 

(17) 



Equation (17) represents the # component of the radiated field. A similar expression 
holds for the 6 component of the field so that once the coefficients b ( * +1) and a ( ^ +1) are deter- 
mined the field may be computed numerically. Additionally, it has been shown that the inte- 
gration over h yields no radiated field except at the stationary points. It should be pointed out, 
however, that the h integration passes through branch points at Uj=0 along the real h axis 
(fig. 2). A plasma sheath differs from a dielectric sheath in that for a plasma the branch 
points are within the region — k <CJie h<Ck , which is the region of integration from which the 
radiation field is obtained. To remove difficulties that arise because of the presence of the 
branch points along the integration path, a small loss term may be added to propagation con- 
stants in the layers, thereby removing the branch points from the real axis. The branch lines 
may then be drawn so as not to cross the real axis and the Fourier transform remains valid. 

k IMAGINARY h 



I I 

I I 

I I 

I ! 

I I 



Tf 7 ?• 



Figure 2 Integration path in the h-plane. 
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3. General Characteristics of the Far Field With Several Plasma Layers 

The general solution for the radiation field when several plasma layers are present is 
difficult to reduce to a useful form because of the size of the expressions and the tedious work 
required to obtain them. As indicated in section 1, the solution for the field when a single 
noncontiguous layer is present requires the reduction of 12 linear equations, and this number 
increases by 4 for every additional layer. Clearly, then, the field expressions are not amenable 
to an extensive quantitative analysis without resort to numerical solutions obtained with the 
computer. Nevertheless, some qualitative and semiquantitative information may be gathered 
about the field pattern without an actual solution of the equations. 

3.1. Polarization of the Field 
The expressions for the radiated field are 

m = — oo J — co 

H,= S r dh\^u,a^mZ ) \Ti,p)\ iM+h * (18) 

7W= _co J-oo L M^ J 

where E^ and H^ represent the electric and magnetic field vectors respectively of the orthogonal 
components of the radiated field. It will be shown in the following section that 

a,n +1) H ( %(u p) is odd in m and 

b { n +1) H ( %(u p) is even in m, 

and since the derivatives of the Hankel functions have the same behavior in m as the Hankel 
functions themselves, eq (18) can be written in the form 

^ J -co 

H,=± e f P dh[-^u a^H^(v p)~] «-«*[«-<"*-**»♦] (19) 

* J-co L M^ J 



where 




(20) 



The field expressions in (19) indicate that the radiated field is in general elliptically polar- 
ized with E<t> given by the upper expression in eq (19) and 2? given by 

Ed =J«LH, (21) 

where H^ is the lower expression in (19). In the planes 0=0 and (t>=w, the E e component is 
seen to be zero and the field is plane polarized. Elsewhere, both components of the field are 
present so that the plasma sheath is seen to introduce a transverse electric field component that 
is not present without the sheath. The order of magnitude of this transverse field, which 
determines the extent of the ellipticity of the radiated field, is difficult to ascertain without 
actual numerical computations, but it is expected that the degree of ellipticity will be small. 

3.2. Proof of the Even and Odd Nature of the Modal Coefficients 

In order to symbolically depict the 4 (iV+1) equations needed to determine the modal 
coefficients, figure 3 is generated, which may be thought of as a representation of the matrix 
that multiplies the line matrix of unknown coefficients to yield the source function P m (h). 
The shaded areas in the illustration are the nonzero sections of the matrix, although particular 
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MODAL COEFFICIENTS , 
(I) (I) (5) (5) 

Orn (|)Arn (|) (2) (3) (4) (5) bm (5) B m 



= 
= P m (h) 



Figuee 3. Symbolical representation of 4(l\-\-l) 
linear boundary matching equations for N = ^. 
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^ H !>4»3> 


' i "f"4 R l Z| ' (u 4'3 1 
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^4 
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Figure 4. Detail of a typical 4-column from, figure 3 showing the multipliers of the modal coefficients for region 4 

at the boundaries p^ and p\. 

coefficients within the shaded section are zero as is indicated in figure 4 which is a detail of the 
shaded region of a typical 4-column of figure 3. 

It may be seen from figures 3 and 4 that the m dependent terms required to compute 
a ( % +1) and 6 ( m +1) are the Bessel and Hankel functions and their first derivatives and the 
quantity m itself. The function P m (h) is independent of m for the line source. In view of this 
the notation below is used 



,(AM-D- 



696-013—63 10 



=/(m, H m , J m ) 

b^=g(m,H m ,J m ). 
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(22) 



First it will be shown that 



then that 



and finally that 



/(m,ff_„,J_„)==(— 1) m j(m, H mj J m ) 

g(7n,H_ m ,J_ m ) = (-l)- m g(m,H m} J m ) (23) 

/(— m, H m , J m ) = —j{m, H m J m ) 

g(—m } H m , J m ) =g(m, H m J m ) (24) 

a { % +1) H^ (u oP ) is odd in m 
6 ( ^ +1) i? ( ^ ) (^ p) is even in m. 

Consider the equations depicted in figure 3. A solution may be found for a { % +1) and 
b { Z +1) by taking the quotient of the appropriate determinants. Every term in the determinant 
in the denominator will contain either a cylinder function or a zero. The determinant in the 
numerator will be missing a column containing cylinder functions since the forcing function will 
replace the appropriate column. The Bessel and Hankel functions and their derivatives behave 
as [Sommerfeld, 1949] 

J- m (x) = (-l) m J m (x) (25) 

so that 

b^ 1) = g(-m,H. m ,J. m )=g(~m,(-irH m X-l) m JnX (26) 



Now, a number common to a column of a determinant may be factored from the determinant 
without changing the value of the determinant. If all the H^^x) and J- m (x) in the expressions 
for a { -ti l) and b^-^ are replaced by the equivalent terms given in eq (25) and then the factors 
(— l) m are canceled from corresponding columns in the determinants in the numerator and 
denominator, the result will be that all such (— l) m terms will cancel but one, left over in the 
denominator. Hence, 



/(m, £U, J. m ) = (- l)-yOn, H m J m ) 
g(m, H_ m , J. m ) = (— 1) - m g(m, H m , J m ) . 



(27) 



MODAL COEFFICIENTS 



(I) (I) (2) 



(5) 
b m 



(3) 



(4) (5) (5) | 
— — ^a m lA m 
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X m x 



X 
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= 

= Pm(h) 



Figuee 5. Symbolical representation of 4(N+f) 
linear boundary matching equations for N = 4 with 
zeros, m, and non-m coefficients placed in evidence. 
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Now, to find the behavior of a ( £ +1) and b { % +1) with m itself, once again consider figure 3. 
If that figure is reconstructed so that zero coefficients, coefficients containing m, and coefficients 
not containing m are placed in evidence, the symbolic representation of figure 5 is obtained for 
the boundary matching equations. Clearly the 4 (tV+1) equations can be reduced by elimi- 
nating one variable at a step starting with a ( ^ ) on the left side of figure 3. Examination of the 
continued application of this process reveals that coefficients that arc even multiples of m are 
always added to even multiples of m and vice versa and, furthermore, that the forcing function 
will always have factors that are of the same even or odd nature as the factors which multiply 
b ( S and B%\ If the process of reducing one variable at a time is continued until the solution 
for 6 ( ;), +1) is found, it is clear that 6 ( ^ +1) has the same even or odd character as P m (h), which 
is even in m, and a ( % +1) has an odd behavior in m. Thus 

/(— m, H m , J m ) = —f(m, H m , J m ) 

g(—m, H m J m )=g(m, H m , J m ) (28) 

or 

/(— m, H- mj J- m ) = —j(m, H-m, J- m ) 

<i(—m, H- m , J- m )=<j(m, H- m J- m ). (29) 

On taking eqs (27) and (29) together, it is found that 

f(-m, H- m J_ m ) = -(-l)-™/(m, H mj J m ) 

<j(-m, II_ m , J. m ) = (-l)- m (j(m, II m , J m ). (30) 

If the expressions in eq (30) arc multiplied by II { ^(u {) p), which has the m behavior indicated 
in eq (25), the desired result is found: 

a^miiuop) = -a^ l) Hir (u oP ) 

V»tPH^oP) =6 < r ,, ^ ) ("oP). (31) 

3.3. Critical Angles of the Field 

The radiation field in the direction d was shown to be contributed by the integration along 
the infinitesimal portion of the real h axis about the point 

h=k cos 0. (32) 

Additionally, it was shown that singularities at the points 

h=kj (33) 

could be removed from the real h axis with the addition of a small loss term to the propagation 
constants in the layers. The contribution to the radiation field of these singular points take 
place in the directions 

cos 0= ±kj/k (34) 

which are exactly the critical angles defined in ray optics at which total reflection occurs. For 
angles of ray incidence beyond the critical angle ray optics predicts no transmission. An 
exact solution to the reflection problem of a plane wave incident at greater than the critical 
angle on an infinite plane plasma slab of finite thickness, however, indicates that power is 
coupled through the slab. The amount of power coupled decreases with increasing angle of 
incidence beyond the critical angle and drops off at a faster rate as the electron density and 
the thickness of the slab are increased. 

For the cylinder problem the solution may be found explicitly for regions defined by 
angles within and beyond the critical angles of the various layers. Within the first critical 
angles, the field considered as a function of 6 depends on the behavior of the various Hankel 
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and Bessel functions of real argument. When the critical angle is reached, in an imaginary 
traversal of the 0=0 plane through angles in 6 starting at 90°, the arguments of some of the 
Hankel and Bessel functions become imaginary. Interpolating from the behavior of plane 
waves propagating through slabs, the radiated fields are non-zero at the critical angles and 
decrease outside the zone defined by these critical angles. In section 4 some computed radiation 
field configurations are shown, but these computations were limited to the region within the 
first critical angles because of the added computational difficulties involved in going beyond 
this noint. Thus an actual determination of the fields in the region beyond the first critical 
angle for the parameter values considered must await further computations. To aid in visual- 
izing the above comments, some typical critical angles are shown in figure 6. 



ZONE OF SMALL 
RADIATION 







CRITICAL ANGLES 
(e>90°) 




Figure 6. Critical angles of the field for a three- 
layer plasma. 



Figure 7. Cylinder encircled by a single plasma 

layer. 



4. Exact Solution for the Single Plasma Layer Problem 



In this section the radiation pattern of an axial magnetic line source on a conducting 
cylinder is considered when a single homogeneous plasma layer which is not contiguous with 
the cylinder is present. This geometrical configuration is of importance because of its obvious 
relation to the problem of propagation through the boundary layer region of the hypersonic 
shock wave. 

From the discussion in section 2 the solution of this problem involves the solution of 12 
linear equations, or actually of 10 nontrivial equations if A ( % and B ( % are excluded in advance. 
The solution for b { %(h) was found, but due to its length only the computed results are presented. 

The expression for the radiation field as given in (17) may be factored into two parts in 
the form 



E*=E°f 



Aperture 



V(z)e ikQZCOae dz 



(35) 



where E% is the field of a point source on the cylinder and the other factor includes the element 
and array factors. A proper terminology for E% is thus "cylinder factor. " This cylinder factor 
was programed for computation on the IBM 7090 for the case of a single plasma layer present 
and for the plasmaless case. 2 The program was established with the electron density and 
cylinder and plasma radii on separate data cards so that these parameters could be changed 
readily. (See fig. 7.) 

As stated above, the solution for b { % is quite lengthy and complicated, as is the machine 
program used to compute the field. In order to develop an independent check on both, a set 
of solutions for the plasma case was taken for decreasing values of electron density and, indeed, 



2 The field with o plasma present is given by [Wait, 1959c J. 
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the field was found to converge to the plasmaless case. The minimum value of electron density 
taken was that corresponding to a relative permittivity of 0.99. The convergence properties 
of the series were also examined through separate print outs and found to be quite good. The 
field generally converged to five figure accuracy in about 14 terms where the outer plasma 
radius, yielding the largest cylinder function arguments, was kept constant at 

67T 

p2= F* 

ft-0 

Computations were made in the plane of the line source and the axis of the cylinder, and 
results are presented in figures 8a through 8g. To insure computational accuracy, fairly 
modest values for the radii were chosen. The cylinder radius was kept constant at 

2tt 
Po=y> 

Kq 

Q ^)7T 47T 

the outer plasma radius as above, and the inner plasma radius was taken at -j— and j-j 
a difference of a quarter free space wavelength. The electron densities correspond to 

-=0.40, 0.80, 0.90, 0.95, and 0.99. 

The chief characteristic of the figures is the extremely sharp variations in field magnitude 
at certain angles. In order not to miss the sharp fluctuations, the field was computed at }i° 

intervals for the - =0.4 and 0.8 cases, and even here the mesh was often not fine enough. In 

figure 8c (- =0.8, pi='-j— > 0=42.5° ) the field rises to a peak 25 db above the free space field 

in a beam only %° wide. This beam and other beams were checked at Ko° intervals, but in 
most cases a simple graphical interpolation was used to complete the curve between ){° in- 
tervals. The presence of these sharp beams implies the existence of leaky waves that travel 
down the cylindrical, lossless, infinite structure and produce highly directive beams which add 
to the other fields. The nature of the radiation field of leaky waves traveling along such a 
structure and the effect which might be expected if the plasma were lossy and of finite length 
are discussed in the next section. 

A second characteristic of the patterns is that the field variations take place around an 
average which is approximately the plasmaless field. The total power radiated in the region 
between the critical angles appears to be about the same whether the plasma is present or not 
and the field near broadside remains fairly constant for the electron densities considered here. 
A low-power antenna with a narrow broadside pattern might be expected to operate fairly well 
in such an environment, but wide and/or scanned beams would be severely distorted. 

5. An Approximation Technique for Quasi-Cylindrical Problems 

It was pointed out in the introduction that the aspect of the plasma sheath problem which 
is of practical interest for mapping and related applications is the one in which the coupling 
between the rf energy and the sheath is fairly light; that is, the case in which the maximum 
relative permittivity in the sheath does not differ by very much from unity. This aspect of 
the problem suggests that the general class of approximation techniques developed from 
integral equation solutions for the field are applicable. 

In the following paragraphs an exact integral expression will be developed for the field 
and the feasibility of approximating the solution to problems related to the cylindrical plasma 
sheath problem will be examined. 
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Figure 8a. Sheath relative permittivity = 0.9 9; sheath 
inner radius = 3. 57r/k . 



1 

SINGL 

PLAS 


E LA 
MALE 


YER f 
SS F 


IELD 
ELO 












1 






































































1 












l^P. 


\~~^ 


















1 ; 


A 


t\ 














































6*^ 


A 


L 


I 




































\ 







90 85 80 75 70 65 60 55 50 45 40 35 30 25 20 15 10 5 



1 1 1 1 

SINGLE LAYER FIELD 

PLASMALESS FIELD 




























| 












































































_i 






4A 




n 

cjl 




ftlfjji i\ \\ 


















r 


v 




















1 




* IT 
1 


r 


If" 




































ii 









90 85 80 75 70 65 60 55 50 45 40 35 30 25 20 15 10 5 

> 



Figure 8b. Sheath relative permittivity =0.95; sheath 
inner radius = 3 .5 t /ko. 



Figure 8c. Sheath relative permittivity = 0.90; sheath 
inner radius = 3. 5 7r/k 






SINGLE LAYER FIELD 
PLASMALESS FIELD 
































\ 


















1 






It 1 . < 










_ — 


r _ J 


<\ 


ii 


A. 1 


jj 


jM 


h 












r 


v 


- 1 ~ 


"^ 








» 












I H 


s 






















1 1 


' 







SINGLE LAYE 

PLASMALESS 


? FIELD 
FIELD 
























1 1 
























jj 


I 


\ 








1 








L_ > 






—A 
— P~" 


At 


-#^4-J^LiA 














[ 














Vrr- 

v 1, 


























1/ 





75 70 65 60 55 50 45 40 35 30 25 20 



75 70 65 60 55 50 45 40 35 30 25 20 



Figure 8d. Sheath relative permittivity = 0.80; sheath 
inner radius = 3 .5ir /ko. 



Figure 8e. Sheath relative permittivity = 0.80; sheath 
inner radius = 4?r/ko. 
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Figure 8f. Sheath relative permittivity = 040; sheath Figure 8g- Sheath relative permittivity =0.40; sheath 

inner radius = 3.5-kJ^. inner radius = ^Tr/ko. 

Figure 8. Radiation fields through a single-layer cylindrical plasma sheath computed in the plane of the axis of 

the cylinder and the line source. 
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5.1. Development of the Approximate Field Expression 
Maxwell's equations in isotropic media, in which e is a continuous function of space, are 
VXE=-jo)fxH (a) V.eE=0 (c) 

VXH=jcoeE (b) V.H=0. (d) (36) 

The curl of (36a) is taken and curl H substituted from (36b). This results in 

VXVXE = «VE (37) 

VXVXE=V(V.E)-V 2 E 
V.eE=eV.E + (V€).E (38) 



When the vector identities 

and 

are employed, (38) reduces to 



where 



(V 2 +F)E = - vi(Ve-E) (39) 



fc 2 = w 2 /*€. (40) 



Now a scalar Green's function is defined; it satisfies the scalar Helmholtz equation with a 
point source inhomogeneity, i.e., 

(V 2 +^)ff(r,r')=S(r-r') (41) 

where r and r' are coordinate vectors. To obtain an integral expression (39) is multiplied by 
G, and (41) by E. These are then subtracted, integrated over space, and the coordinate 
variables are reversed to obtain 

E(r)= f(EV*G-GV 2 E)dv'- fc7vi (Ve . E)dv'+ f(k 2 -k 2 )GEdv'. (42) 

The first term in (42) may be further reduced to a surface integral by applying the scalar 
Green's theorem. Additionally, if a Green's function is selected so that it vanishes on the 
cylinder surface, eq (42) becomes 

E(t)= fE(V#. ds)+ f(k 2 -k 2 )GEdv'~ fc7V^ (V*. E)dv\ (43) 

Equation (43) is an integral equation for the electric field which has little utility when 
an exact solution for the field is required. Even for the case in which e=e Q throughout the 
region external to the conducting cylinder, i.e., when the second and third terms in (43) are 
zero, a direct solution cannot readily be found because the surface integral in (43) must be 
taken over the radial component of the field on the cylinder as well as the source field in the 
aperture. Nevertheless, a good approximation to the radiated field may be found by taking 
the surface integral over the aperture alone [Silver, 1949]. This approximation is especially 
true in regions near broadside. 

For the inhomogeneous case the effect of the plasma on the radiation field within the 
critical angles is a perturbation on the fields when no plasma is present. Further, the first 
term in (43) represents the plasmaless field. The error in this assumption is once again neglect 
of the radial field at the cylinder, which is altered by the presence of the plasma. The second 
and third terms in (43) thus are taken to represent the field perturbation due to the plasma 
and, in fact, the third term will be neglected; it is essentially an integration over the localized 
charge density in the plasma. That is, the third term is considered a less significant term than 
the second which represents an integration over the currents in the plasma. 
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The above comments are represented by (44) and (45) 

-■^total ~*Vo plasma ~T -^plasma v*4/ 

E vl&sma ~f(k 2 -k 2 )GEdv' (45) 

5.2. Integration Over the Plasma 

The scalar Green's function defined by (41) and the condition that (7=0 at the cylinder is 
[Wait, 1959c] 

G=^ jt m f^ dhH m (u oP ) [j.CnoPO-Hff^PO h^Tpo) ] e- i[m{ *~*' )+h(z - z ' )] (46) 

where the primed coordinates refer to the space containing the plasma and the unprimed co- 
ordinates refer to the field space. To find any particular far field component the function 
(k%—k 2 )0 is multiplied by the components of E in the plasma parallel to the far field vector. 
For example, to find the </> component of electric field radiated by the plasma 

s W.= f «-F)(7E(r')-e>' (47) 

U plasma 

where e^ is the unit vector parallel to the far field vector. Of special interest is the plane 0=0. 
In this plane, (47) and hence the field radiated by the currents in the plasma reduces to 

s ** m =f (kl-Jc*)G[co S <t,'E r +sm<t>'E p ,W. (48) 

%J plasma 

Now, the term inside the square brackets in (48) may, in general, be expanded in a Fourier 
series in <j>' and a Fourier integral in z' to yield 



cos 



$'E r +sm<j>'E p ,]=Jt r dla n (l,p')e- i{n *' +lz '> (49) 

— oo J —co 



where n and / are the transform variables. This form for the field is of course the same form as 
was employed previously for the exact solution to the cylinder problem. Substitution of (49) 
and (46) into (48) results in 

dz' P 'd P ' d4> Z! £ dh\ dl 

-oo J PO JO 77l= — co n = — oo J - co J —co 

■f (K-t?)e- ,l -~* ,+hl *-" )] H%(u oP ) [j m (u oP ')-H%(u oP ') i$* Po \ ~\-c t ,(l, p ')«-«-+'+uo |- (50) 

A cylindrically symmetrical plasma sheath has been assumed so that the integration on <$> r 
may be carried out directly to yield 

#w= 2 *- r ^ r p,(i >' s r dh r * 

J -co J PQ m= — co J - co J —co 

•/(«-*•)«-*<— '>a g>(«op) [j m (u p')-H™(uy) h%l^t) \ am{1, p ' )e ~ ilz '}- (51) 

In a similar manner the integration on z r may be carried out, but here there is the possibility 
that (kl— k 2 ) varies with z'. In that event the z' integration yields the Fourier transform of 
(kl—lc 2 ) in the transform variable (h—l) . The integration over I will then yield the convolution 
of the transform of (kl— k 2 ) witha m (Z, p'), which is denoted by p m (h, p f ). Then (51) reduces to 
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JPO 771= — co J -oo 

■/a.(a, p')e- iftj ^'(woP) [j„(w„p')-tf;» 2 w') h!^) ]}- (52) 

The integration on A may be carried out for large p as in (17) to yield 

^*-=-^ £ ^f"p'^' s -ffl.(A,p / )«*-' A -rJ r -(tw»')-fl?(«op / )^^i} 

(53) 

5.3. Characteristics of the Radiation Fields of the Plasma Currents 

Application of (53) to the approximate solution of the type of cylindrical and quasi- 
cylindrical problems of concern here is by no means simple. The prime source of difficulty 
lies in the extremely limited class of functions @ m (h, p) [Sonine, 1880] which, when a factor in 



the integrand 



p'pjh, p') \j m {u nP ')-ii^{uy) ^riy°\ 1> 

L "m \ u qPo)J 



yield an integral which is a known tabidated function. Indeed, the cylinder functions which 
are solutions to the exact boundary value problem when taken as factors of p m (h, p) do not all 
yield known functions for the p f integration in (53). Considerable investigation of the p' 
integration has been made. The net result of this investigation is that the integration over 
p yields a slowly varying function in 6 that represents a small addition to the free space field. 
It was expected that this small perturbation would be the extent of the pattern degradation 
within the critical angles, but the fields computed for the single layer case indicate more sig- 
nificant perturbations on the field configuration than can be attributed to the terms which 
result from the p integration. In general, the p integration will yield more significant pertur- 
bations as the thickness of the sheath increases. This fact may be gathered intuitively by 
considering the equivalent antenna aperture size represented by the sheath thickness and an 
image current sheath in the cylinder. 

The perturbations which are evident in figure 8 are attributable to the presence of leaky 
waves that propagate along the cylinder and sheath and give rise to highly directive radiation 
patterns. A leaky wave excited at z' = z u and traveling outward along the C3 T linder with 
propagation constant v } and attenuation constant yj in the form 

[cos 0'£V+sin 4>'E,.]=f(p')e ti 'i +i W- t "i (54) 

will radiate a directive field with the beam pointing direction given by 

cos 0= ± v f (55) 

and with the beam width determined by ~ as will be indicated in the following paragraph. 

To compute the field radiated by the leaky waves, consider the beta function for such 
waves when the sheath is uniform in the z' direction 



r (I /w^- /w J (p/) / l I l \ 

where the radial dependence of the function 



(56) 



[cosfEr+siiKp'E,,,] 
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for the leaky waves is included in the f unction f m .(p') an d the index j is subscripted to the m 
to indicate the various leaky waves in each mode. On substituting (56) into (53) there results 



e -iJcR ( \ 1 \ 

'{I" p ' dp '^ (p ' )elm;f [^(^pO-flffC^ 



' HX?(u oPo ) 



jfh = kQ<. 



(57) 



The integral within the second braces is generally a slowly varying function of 6 as was dis- 
cussed above. The integration over z" which is implicit in the second bracket in (57) and 
yields the field pattern of a linear array may negate the above statement, but only the cylinder 
factor, which is the field due to a point source, is considered here. 

The first bracket in (57) yields the rapidly var}dng dependence of the field on 0. The 
first term reaches peak magnitude at 



and the second term at 



h=Jc cos d=-~Vj 
h=k cos 6=Vj 



(58) 
(59) 



which are at symmetrical points about broadside. It will be recalled that the field in (57) is 
to be linearly added to the free space field in this approximation so that the real and imaginary 
parts of the rapidly varying functions in (57) are of interest. Figure 9 is a sketch of the 
real and im aginary parts of 

-1 



-jj+iih—uj)] 



near 



=cos 






(60) 



Examination of figure 8 will indicate that the perturbations on the free space field can be 
accounted for by the complex addition of functions of the type shown in figure 9 to the free 
space field. 

The magnitude of the leaky waves may be approximated with reasonable accuracy from 
figures 8a through 8g by considering that the peak or null of each sharp field perturbation is 
the beam pointing direction of the field radiated by the jth leaky wave. The propagation 
constants, v j} and the attenuation constants, y j} can also be approximated using the peaks 
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Figure 9. Sketch of the real and imaginary part of 
a field radiated by leaky waves for small angles near 
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and, in addition, the widths of the perturbations. The analytical determination of vj and jj 
necessitates the solution of an extremely complicated transcendental equation which involves 
the cylinder functions. This latter method has not been carried out, but it is suggested that 
the approximations represent a useful method by which these eigenvalues can be obtained. 

Without an analytical expression for the propagation and attenuation constant of the 
leaky waves, predictions of the effect on the radiation pattern of non-zero losses in the plasma 
and a finite length of plasma are difficult. One might expect, however, that non-zero losses 
in the plasma would cause an increase in the attenuation factors jj which (57) indicates would 
mean a broadening and lowering of the perturbations on the free space field. Similarly, 
limiting the sheath to finite length would alter the p m (h, p') function in a form equivalent to 
the decrease in aperture size in the z direction and introduce a source field at the disconti- 
nuities. Both of these factors would tend to smooth out the field perturbations. The field 
pattern in the presence of a hypersonic plasma sheath would thus not be expected to have the 
sharp perturbations of this simplified model. 

6. Conclusions 

It has been demonstrated that quite a significant distortion of the field pattern is effected 
by the presence of a lossless plasma sheath. The chief characteristics of this distortion are: 

1. The presence of critical angles in the field which reduces the power radiated in certain 
directions so that a blackout zone is introduced in which the field probably is extremely small. 

2. The presence of narrow but large peaks and valleys in the field between the critical 
angles which are attributable to the generation of leaky waves that travel along the cylindrical 
structure. 

3. The generation of a transverse field that changes the plane polarized radiation field 
when no plasma is present to an elliptically polarized field in every direction except the planes 
</>=0 and <j)=ir. 

4. The presence of slowly varying perturbations on the field which are on the order of 
several db due to the non-leaky wave fields in the plasma. 

The effect of these distortions on the total field pattern depends on the aperture distribu- 
tion. In (35) this dependency is made evident: the total field is indicated to be the product 
of the element, array, and cylinder factors. Only the cylinder factor is shown in the curves. 
For array factors which yield narrow, broadside beams with very low sidelobes, the distortion 
under the parameter values discussed in this paper will be represented by increase or decrease 
of only several db in pattern gain of the antenna. When the beam is scanned be}^ond a few 
degrees, however, or when the array factor is a broad beam, a great deal of distortion occurs 
for even modest electron densitites in the sheath. 
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